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ABSTRACT 

We present a study of the flickering activity in two nova like systems KRAur and 
UUAqr. We applied a statistical model of flickering simulations in accretion discs 
based on turbulent angular momentum transport between two adjacent rings with an 
exponential distribution of the turbulence dimension scale. The model is based on a 
steady state disc model which is satisfied in the case of hot ionized discs of nova like 
cataclysmic variables. Our model successfully fits the observed power density spectrum 
of KRAur with the disc parameter a = 0.10 — 0.40 and an inner disc truncation radius 
in the range Rin = 0.88 — 1.67 x 10^ cm. The exact values depend on the mass transfer 
rate in the sense that a decreases and i?in increases with mass transfer rate. In any case, 
the inner disc radius found for KRAur is considerably smaller than in quiescent dwarf 
novae, as predicted by the disc instability model. On the other hand, our simulations 
fail to reproduce the power density spectrum of UU Aqr. A tantalizing explanation 
involves the possible presence of spiral waves, which are expected in UU Aqr, because 
of its low mass ratio, but not in KRAur. In general our model predicts the observed 
concentration of flickering in the central disc. We explain this by the radial dependence 
of the angular momentum gradient. 
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1 INTRODUCTION 

Cataclysmic variables are interacting binaries with a late 
type secondary star overfilling its Roche lobe and transfer- 
ring mass from the inner Lagrangian point Li toward a white 
dwarf primary. When the white dwarf in not magnetic (or 
weakly magnetic i.e. an intermediate polar) an accretion disc 
is formed. Disc matter is transported inwards whereas angu- 
lar momentum is transported outwards. Accretion is a com- 
mon phenomenon of a broad family of astrophysical objects. 
X-ray binaries are the most similar accreting systems to cat- 
aclysmic variables. They possess a black hole or a neutron 
star as central objects. Supermassive black holes are also 
the central engine which powers active galactic nuclei. Ac- 
cretion of matter through a disc is the source of brightness 
variations on a wide range of time scales and energies. The 
time scale of the variability scales with the mass of the com- 
pact object and this is the reason why cataclysmic variables 
provide excellent laboratories to study accretion physics. 
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The most characteristic signature of accretion is the 
flickering activity. These brightness variations, with ampli- 
tudes in the range of 0.1 — 1 magnitude, last from seconds to 
tens of minutes. There are four basic physical interpretations 
for the flickering activity that we know of; 1) interaction of 
the gas stream from the secondary with the outer edge of 
the disc , 2) turbulent flow in the accretion disc, 3) mag- 
netic dissipation/reconnection in the disc corona and 4) in- 
teraction of the inner edge of the disc with the central white 
dwarf, the so-called boundary layer. Furthermore, a c e llular - 
automaton model was proposed by lYonehara et al.l (|l997l ) 
to describe the timing properties of flickering in cataclysmic 
variables. In this model light fluctuations are produced by 
occasional flare like events and subsequent avalanche flow 
in the accretion disc atmospheres. Flares are assumed to 
be triggered when the mass density exceeds a critical limit. 
There have been several attempts to localize the source of 
flickering. For example, enhanced flickering activity during 
th e orbital hump before the eclipse of U Gem was observed 
by IWarner fc Nathed l|l97ll ). This allowed them to identify 
the flickering source with the interaction of the gas stream 
with the outer edge of the disc (hot spot). On the other 
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hand, iBruchI l|l992l ) studied energies and variability time 
scales in a large sample of cataclysmic variable and sym- 
biotic systems and concluded that both the disc and the 
boundary layer are responsible for the flickering. The com- 
bination of the inner disc and hot s pot is the sou rce of flicker- 
ing variability in the case of Z Cha lBru ch'1996'') and HT Cas, 
V2051 Oph, UXUMa, IP Peg (|Bruchl2 000'). Eclipse mapping 
of V2051 Oph identified the inner disc as the source of high 
frequency flickering and the stream overflow as the sourc e 
of low frequency variability (|Baptista &: Bortolettd |2004| ) . 
The inner disc is also r esponsible for the flicker ing in the 
symbiotic system T CrB ijZamanov fc Bruchll998l l . This was 
studied and conflrmed by statistical model ing of turbulent 
accretion disc flow by iDobrotka et al.l (|201Q ) . The authors 
explained the central concentration of flickering by the large 
number of turbulent eddies required to transfer angular mo- 
mentum from the inner disc radius outwards. Therefore, the 
accretion disc appears as the most common source of flick- 
ering in cataclysmic variables and related object such as 
symbiotic systems, specially its inner parts and boundary 
layer. 

Nova lik es are a subclass of cataclysmic variables (see 
IWarneilll995l for review) with high mass transfer rates. Their 
discs are hot and fully ionized most of the time which means 
that the mass accretion rate is constant through the whole 
disc (i.e. steady state). This is opposite to dwarf novae where 
the mass transfer rate is lower, the mass accretion rate is 
not constant and intersects the critical value of thermal sta- 
bility somewhere in the disc at some time. This drives a 
viscous-thermal i nstability which triggers the d warf novae 
phenomenon fsee lKato et al.iri998l . lLasotall200ll for review). 
Dwarf novae discs are not in steady-state and switch be- 
tween hot (ionized) and cold (non ionized) states with the 
former presenting smaller disc truncation radii. A small in- 
ner disc truncation radius is expected also in the case of nova 
like systems. The inner disc truncation radius can be pro- 
duced either by weak magnetic fleld (i.e. not strong enough 
to qualify the system as an intermediate polar) or some kind 
of evaporation (see e.g. iMever fc M cvor-Hofmcistcr 1994). 
Th e case of ionized steady sta te discs is well described by 
the lShakura fc Sunvae^^ (| 19731 ) accretion disc model. 

Extensive work has been devoted to locali ze the source 
of flickering in cataclysmic variables (see e.g. BruchI 19961. 
Zamanov fc Bruchll 19981 . iBruchlbOOd iBaptista fc Bortolettd 



2004 ). but we know little about the physical phenomenon 



responsible for its variability. The location of flickering in 
the inner disc and boundary layer incorporates three pos- 
sible physical processes: turbulence, magnetic reconnection 
and interaction of the inner disc with the central star. In 
this paper we are trying to test whether turbulence is the 
main mechanism responsible for flickering by applying the 
statistical model of t urbulent accretion flow developed in 
iDobrotka et al.l (|201(]| ). The model is applied to the flickering 
activity observed in the nova likes KRAur and UUAqr. In 
Sec. [2] we summarize the observational flickering properties 
of KRAur and UUAqr. In Sec. [3] we retrace the statistical 
turbulence model of Dobrotka et al. (2010). The uncertainty 
of the measured parameters is discussed in Sec. 3] The fl- 
nal results are presented in Sec. O discussed in Sec. [6] and 
summarized in Sec. [T] 



2 OBSERVATIONAL RESULTS 
2.1 KRAur 

KRAur is a nova like sys tem of VYScl ty pe. Its flicker- 
ing activity was studied bv lKato et al.l (|2002l ) during a high 
state. The authors used 17 observations, with a mean cover- 
age of 9000 s and 50 s time resolution, to calculate the power 
density spectra (PDS). The PDS of the flickering is well de- 
scribed by a broken power law, characteristic of red noise. 
The power law has a slope with index of -1.63 and a cut off 
frequency at 31.6 d~^ (or 1.5 d~^ in logarithmic units). 

Quasi-periodic variability was detected on some nights 
with a characteristic time scale of 10 — 15 min. No super- 
humps or other coherent signals were observed. 



2.2 UUAqr 

Photomet ric light curves of the n ova li ke UUAqr were ob- 
tained by IBaptista fc Bortolettd (j2008') with 10 s time res- 
olution and a mean duration of 8000 s per run. The PDS 
shows a power law with index -1.5 and a break frequency 
at ~ 12.96 d"'^ (or 1.11 in logarithmic units). The eclipse 
mapping technique was applied to analyze the short term 
variability. Low and high frequency flickering maps are dom- 
inated by emission from two asymmetric arcs simi lar to the 
case of the dwarf nova IP Peg l|Baptista et al.ll2002l ). The in- 
terpretation of the flickering activity in UU Aqr is then sim- 
ilar to IP Peg, i.e. tidally-induced spiral shock waves in the 
outer regions of a large accretion disc. Baptista et al. (2002) 
have suggested that the flickering in UUAqr is caused by 
turbulence generated by the disc-spiral wave interaction. 



FLICKERING MODELING AND 
SIMULATIONS 



IDobrotka et al.l (|2010l ) developed a statistical method to 
model the flickering light curves using turbulence. The 
method is based on angular momentum transport between 
two adjacent concentric rings in a disc. The angular momen- 
tum is transported by discrete turbulent spherical bodies 
with a dimension scale x described by an exponential dis- 
tribution function. The dimension scale is proportional to 
the radial propagation of the turbulent body in the accre- 
tion disc. The maximum dimension is the scale height of the 
disc, while the minimum is limited by zero. 

The basic idea is that the sum of the angular momen- 
tum difference AZ(r) (where r is the distance from the disc 
center) of all spherical bodies in two adjacent rings must be 
equal to the global angular momentum difference between 
the two rings AI/(r), i.e. Ai(r) ~ AI/(r). Since the scale 
height, density, radial and tangential velocity of the matter 
in the rings is a function of r, the angular momentum dif- 
ferences must be correlated with a certain parameter k{r) 
which depends on the radial position of every pair of rings, 
and is given by 



fc(r) 



AL(?-) 
AZ(r) ■ 



(1) 



Therefore, k{r) is a measure of the quantity of events be- 
tween two adjacent rings at radial distance r. This param- 
eter is needed to estimate the number of flares with dura- 
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tion t ~ x/vi{r), where Vr{r) is the radial viscous velocity. 
The number of events with a certain duration is summed up 
through all pair of rings at distances n 



CW = Vfc(n)/(n,t). 



(2) 



where /(ri,t) is the exponential distribution of turbulent 
eddy sizes with t ~ x/vr{r) at distance ri from the disc 
center and n is the number of ring pairs. By selecting the 
number of flares with durations t between a minimal and 
maximal value and a sampling time step, we get the final 
histogram used in th e synthetic light curve generation (see 
iDobrotka et al.ll2010l for details). 

The disc ra d ial profile is based on the 
IShakura fc SunvaevI l| 19731 ) model. Therefore, the disc 
must be in the hot ionized steady state , characteristic of 
dwarf novae outbursts (see iLasotal 1200 ll for review) . Nova 
like systems are adequate for such modeling because of 
their high mass transfer rates which tend to keep them in 
the hot branch of dwarf novae activity (see I WarneJI 19951 for 
review) . 

We used this method to test whether the turbulence 
scenario is able to fit the observed PDS with realistic values 
of the disc a parameter and the inner disc truncation radius. 
The outer disc radius flo ut is taken as halS radius of the 
primary Roche lobe Rli l|Paczvriskilll97jl : 



Rli = 0.462 a 



Ml 



1/3 



(cm). 



(3) 



. Ml + M2 

where Mi and M2 are the primary and secondary mass re- 
spectively and a is the distance between the two stars cal 
culated from Kepler's Third law 

Mi-^'^' 



a = 3.53 X 10' 



(4) 



q = M2/M1 is the mass ratio and Porb the orbital period. 

The critical mass tran sfer rate for a disc to be in a hot 
stage is given by equation (jHameurv et al.lll998l ) 



Mtr = 9.5 X lO'^a"-''' 



Ml 

Mo 



U0i0( 



Using a values from 0.01 to 1.00 we obtain a critical mass 
transfer rates of 7.8 - 8.2 x lO^'* gs~^ for both KR Aur and 
UU Aqr. By increasing the outer disc radius to ~ 0.9 x _Rli 
the critical mass transfer rate rises to 3.8 — 4.0 x lO^^gs"^ 
for KRAur and 3.8 - 3.9 x 10^^ gs"^ for UUAqr. Therefore, 
we decided to take 8 x lO^^gs"^ as a lower limit to the 
mass transfer rate for a steady disc. Simulations were also 
computed for 1 x 10^'', 5 x lO^'' and 1 x 10^* gs'^ The 
model input parameters and the PDS values, as reported in 
the literature, are summarized in Table [T] 

As mentioned above the dimensionless parameter a and 
the inner disc radius Rin are the free parameters of our simu- 
lation. We calculated PDS power law indexes {pi hereafter) 
and cut off frequencies (/ctf hereafter) for the parametric 



^ The tidal perturbation of the secondary truncates the outer 
disc radius at ~ 0.9 X i?Li and hence provides an upper limit to 
Rout, but our simula tions are not sensitive to this parameter (see 
iDobrotka et al.ll201(]|) 



Table 1. Measured and calculated parameters. Mi is the primary 
mass in solar units, q is the mass ratio, pi index is the power 
law index as measured from the PDS, /ctf the cut off frequency 
in mHz, Porb the orbi tal period in hou rs, the white dwarf 

radius estimated from iNauenberd l ll972l ) in 10® cm and -Rout the 
outer disc radius as half of the primary Roche lobe in units of 
10"' cm. 



Object 


Ml 


g 


pi index 


log(/ctf) 




(Mq) 






(d-i) 


KRAur 


0.59^ 


0.6^ 


-1.63^ 


1.50^ 


UUAqr 


0.67^ 


0.33 


-1.50" 


1.11* 


Object 


^orb 




^Jout 






(h) 


(10® cm) 


(IQiO cm) 




KRAur 


3.907^ 


0.87 


1.87 




UUAqr 


3.900^ 


0.80 


1.95 





Ref erences: 

1 - IShafted lll983ll 



. 2 - iKato et a.1.1 iWQ'j) . 3 - iBaptista et al.l 

1I1994I) . 4 - iBaptista fc Bortolettd | |200S|) 



space a — Rin wi th a resolutio n of 20 x 20 using the Lomb- 
bcare le method (|Scarglelll983 ). The a parameter is allowed 
to vary from 0.01 up t o 1.00 which are typical valu es for cata- 
clysmic variables (see IWarneill 19951 . lLasotall200ll for review), 
whereas Rin varies from the vicinity of the white dwarf ra- 
dius (i.e. 1 percent larger) up to 2.0 x 10^ cm in 20 steps. 
T he white dw a rf rad ius is estimated from the formula given 
m iNauenbergj (|l972| ). For every combination of parameters 
we calculated a synthetic light curve and its PDS. The syn- 
thetic light curves were sampled in identical ma nner as the 
observ ation s of KR Aur and UU A q r repo rted in lKato et al.l 
(|2002| ) and IBaptista fc Bortolettol l|2008l ). These define the 
frequency limits and the resolution of the synthetic PDS. 
We repeated the simulation 1000 times and then calculated 
the mean PDS. The resulting mean PDS is binned into 20 
equally spaced bins and fitted with two power laws. The 
break frequency f^ti and the red noise slope pi are the main 
output parameters. 



4 UNCERTAINTY ESTIMATE 

Unfortunately we have no information about the measure- 
ment errors of the observed PDS parameters used in this 
work. The uncertainty must be taken into account because 
of the li mited number of d ata used in the mean PDS cal- 
culation. Kato et al. [toot only used 17 observations while 
IBaptista fc Bortolettol (|200^ ) used 31. This is considerably 
lower than the 1000 synthetic runs used in our simulations. 
Furthermore, the length of the observed light curves is vari- 
able and some are very short. Our simulated data sets have a 
constant duration in order to get confident PDS parameters. 
This can have an impact in the observed PDS properties and 
must be taken into account in the analysis of the deviation 
between the observed and simulated data. 

In order to estimate this, we performed simulations of 
17 synthetic light curves for KR Aur with identical duration 
and sampling as described in Sec. l2.1l using an ad hoc model 
with Mtr = l.OxlO^'^gs'S a = 0.20 and Rin = 1.0x10® cm. 
A mean PDS was calculated from the synthetic light curves 
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-1.59 -1.58 -1.57 



-1.9 -1.3 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 

power law index 

UU Aqr - simulations 




-1 .975 -1 .97 -1 .965 -1 .96 -1 .955 -1 .91 



-1.95 
/ index 



KR Aur - simulations 

' 2500 
2000 
1500 
1000 
500 



1.55 1.57 1.59 1.61 1.63 



log(cut off frequency) [1/d] 

UU Aqr - simulations 




1.81 1.82 1.33 1.84 1.85 1.3e 



lDg(cut off frequency) [1/d] 



Figure 1. Histograms of tiie simulated PDS parameters obtained using an ad hoc model with Mtr = 1-0 X lO^^gs"^, a = 0.20 and 
iJin = 1.0 X 10^ cm. Left panels show the histograms of the power law index for KR Aur (top) and UU Aqr (bottom) whereas right panels 
present the histograms of the cut off frequency in logarithm units (also KRAur on top and UU Aqr at the bottom). The dotted lines show 
the Gaussian fits to the histograms. Seventeen synthetic runs were produced to generate the PDS of KRAur and 31 runs for UUAqr. 
The histograms in the inset panels were produced using 1000 synthetic runs. The simulations were repeated 10000 times in each case. 



and the process was repeated 10000 times. The resulting 
histograms of the power law index and the logarithm of the 
cut off frequency are depicted on the upper two panels of 
Fig. [TJ together with their best Gaussian fits. The simulated 
PDS parameters cluster at the mean values pi = —1.59 and 
log(,fctf) = 1.58 d~^ with 1-(T of 0.06 and 0.12 d"^ respec- 
tively. We take this l-cr values as our estimate of the error 
in the measured PDS parameters. For comparison, we ob- 
tain pi = -1.580 ± 0.004 and log(/ctf) = 1.60 ± 0.01 d'^ 
when increasing the number of synthetic light curves to 
1000. The corresponding histograms are shown in the in- 
set panels to Fig. [T] A similar analysis was also performed 
for the case of UUAqr using 31 syn thetic light curves per 
mean PDS as in the observations of iBaptista fc Bortolettd 
l|2008l ). The duration and resolution of the synthetic light 
curves is identical to the mean values quoted in Sec. 12.21 
The PDS parameters cluster at the mean values pi = —1.97 
and log(/ctf) = 1.85d~^ with l-cr widths of 0.04 and 0.09 d"^ 
respectively. For comparison 1000 synthetic light curves per 
mean PDS yields pi = -1.961 and log(/ctf) = 1.836 d~^ 
with l-cr widths of 0.004 and 0.007 d"^ respectively. The his- 
tograms of the PDS parameters for UUAqr and their best 
Gaussian fits are depicted in the lower panels of Fig. [T] 



5 RESULTS 

Fig. [2] presents the simulated PDS parameters in the para- 
metric space a — Rin using a grey scale. Only the case for 
Mtr = 8 X lO^^gs"^ is shown for illustration. The isocon- 
tours indicate the measured values (middle line) and the 
±l-cr errors as estimated in Sec. 3] No contours are visible 
in the case of UUAqr because they lie outside the plotted 
a parameter range 0.01 — 1.00. In this paper we are mainly 
interested in the a range 0.1 — 0.5. The lower limit is set by 
the minimum a parameter required by the disc instability 
model (DIM) for hot discs while the upper limit is some- 
how arbitrary. Values higher than 0.5 are theoretically pos- 
sibl e but previous works s uggest that they are unlikely (see 
e^ 



; put previous worKs s uggest tnat tney are unuKely ( see 
Schreiber et"al]|2003l . ISchreiber et"al]|2004 iKing et al.l 



l2007i r Simulations of hot discs in cataclysmic variables typ 
ically yield a = 0.1 — 0.2 whereas observations constrain a 
to the range 0.1 — 0.4. 



5.1 KRAur 

The first clear result of our simulations is that the power law 
indexes measured from the PDS are in the a interval 0.1 — 1.0 
as required by hot discs of nova like systems in the high state. 
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KR Aur - power law 



KR Aur - log(cut off frequency) 




log(a) 

UU Aqr - power law 



-1.5 -1 -0.5 

log(a) 

UU Aqr - log(cut off frequency) 




log(a) 



SxlO^^gf 



Figure 2. PDS parameters for KR Aur (upper panels) and UU Aqr (lower panels) in the parametric space a — Ri^ for Mt, 
Left panels represent the fitted power law and the right panels the logarithm of the cut-off frequency Contour lines indicate the mean 
value of the parameters and their itl-cr limits. For UUAqr these lie outside the calculated range a = 0.01 — 1.00. 



Table 2 . Extreme a and values allowed by the overlap of the 
ztlcr contours in Fig. 3. Different mass transfer rates are consid- 
ered. 



Mtr 


min a 


min J?in 


max a 


max i?in 


(lOi^gs-1) 




(10^ cm) 




(10" cm) 


0.800 


0.10" 


O.SS*" 


0.40 


1.56 


1.000 


0.10" 


O.SS*" 


0.40 


1.60 


5.000 


0.10" 


0.88'' 


0.33 


1.64 


10.00 


0.10" 


0.88'' 


0.25 


1.67 



" Minimum a parameter allowed by theory. 

'' Minimum inner disc radius as set by the radius of the white 
dwarf. 



On the other hand, the cut off frequency points to low a 
values (see upper right panel in Fig. [2]). Fig. [2] also shows 
that the simulated power law index allows untruncated disc 
radii while the cut off frequency do suggests a truncated disc 
(i.e. -Rin > 2.1 X 10" cm for a > 0.1). We then decided to 
search for Rin values in the range a = 0.1 — 0.5 for which 
the power law index and cut off frequency overlap within the 
estimated errors and the results are hsted in Table (2] Fig. [3] 
shows the dependence of the PDS parameters with mass 
transfer rate in the i?in — a plane. For the sake of clarity, 
Fig. [4] also shows the evolution of the mean PDS parameters 
without the error intervals. The displacement of the isolines 
towards lower a. values with increasing mass transfer rate is 
clear. 



5.2 UUAqr 

The results of our simulations for the case of UUAqr are 
totally different than for the previous system. For illustra- 
tion we show the case of the lower limit mass transfer rate 
(8 X 10^^ gs^^) in the lower panels of Fig. [21 The observed 
cut off frequency log(/ctf) = 1.11 is obtained for very low 
values of the a parameter beyond the upper left corner of 
the figure (where the simulated value is 1.2, i.e. minimal). 
Therefore, the whole interval lies outside the calculated a 
range. Assuming the same trend of iiin versus a as in the 
KRAur case, very large inner disc radii would be required 
to get a physically realistic a parameter. Therefore, our sim- 
ulations of the cut off frequency do not yield a satisfactory 
result. Similarly, the observed power law index can only be 
reproduced by our simulations for unplausible values a > 1. 
And we note that these conclusions do not depend on the 
assumed mass transfer rate. 

The statistical distribution of the simulated cut off 
frequencies shown in Fig.l (main figures) are different for 
KRAur and UUAqr. In the former case the histograms are 
well fitted with a Gauss function while in the latter case 
the Gaussian fit only gives a rough approximation. Conse- 
quently, our Gaussian-based uncertainties may not be very 
realistic. However, we note that a minimum error of S-a is 
needed to bring the cut off frequency into required (plausi- 
ble) a interval. On the other hand, the power law histogram 
clearly follows a Gaussian distribution and even a lO-cr error 
is not able to bring the power law interval into acceptable 
a values in the UUAqr case. Therefore, we conclude that 



6 A. Dobrotka, S. Mineshige and J. Casares 




our results are robust against a reasonable error amplifica- 
tion and hence no acceptable fit can be obtained for UU Aqr 
with our simulations. 



6 DISCUSSION 

6.1 Model degeneracies 

By fitting our turbulence model to the flickering data we 
find that the best fitted a parameter decreases while i?in 
increases with the selected value of Mtr- How can this be 
explained? The v elocity of matter along the disc radius is 
given by (see e.g. lFrank et al.lll992l '): 



^4/5^3/10^- 



■1/4^-1/4^-14/5 



where 



/ 



r 



1/2 



1/4 



(6) 



(7) 



Equation[6]shows that the radial viscous velocity increases 
(and hence the flaring time scale decreases) with both A/tr 
and a. Therefore, the simulations performed with higher A/tr 
result in a higher v-c and hence a shorter flaring time scale. 
As a result, the simulated statistics of the flickering changes 
and so do the PDS parameters, such as the power law slope 
and break frequency. In order to reproduce the observed 



power law properties, the drop in flickering time scale or 
event duration must be compensated by other parameters 
such as a. In other words, the impact of increasing Aftr can 
be compensated by decreasing the a parameter in order to 
keep the observed power law unchanged. This is the reason 
behind the A/tr — ct correlation. A similar explanation holds 
for the Aftr — Riu behavior. If we increase A/tr the charac- 
teristic time scale of the flickering events will drop and the 
modeled flickering statistics will contain a larger number of 
fast events. In order to eliminate the excess of fast events 
we need to cut down the inner disc, where short time scale 
events dominate. Therefore, we can compensate for the in- 
crease in Mtr by increasing R j^. A correlation betw een a 
and -Rin was already found bv iDobrotka et al] (|2010l '). The 
explanation is again based on the dependence of the radial 
velocity on a. By increasing a the radial velocity also in- 
creases which causes a decrease in the event duration. The 
quantity of short events is then enhanced and the inner disc 
radius must rise to compensate for it (see Dobrotka et al. 
2010 for details). 



6.2 Results compared with DIM predictions 

The study of flicke ring using the statistical model of 
iDobrotka et al.l l|2010l ) yield two different results for the nova 
like cataclysmic variables KRAur and UU Aqr. A successful 
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fit to the observed PDS properties is obtained for KRAur 
but not UU Aqr. In the former case the observed power law 
index and cut off frequency can be reproduced with a trun- 
cate d disc and pa rameter a ~ 0.1, as required by the DIM 
(see lLasotall20oil for review). Accretion disc theory sets a 
lower limit to a ~ 0.1 while the upper limit depends on the 
mass transfer rate used in the model. The higher the mass 
transfer rate, the lower the a parameter. Maximal values of 
the a parameter 0.25 — 0.40 are obtained for mass transfer 
rates in the range 0.08 — 1 x 10^* gs~^. These a values are 
in excellent agreement with the ones frequently used or de- 
rived in simulations ( 0.1 — 0.2 , see e.g. ISchreiber et"al]|2003l . 
ISchreiber et al.|[200'll ) and with the obser ved values . 1 — 0. 4 
of fully ionized geometrically thin discs (|King et al.l 120071 ). 
On the other hand, the inner disc radius is allowed to vary 
from nearly the white dwarf radius up to 1.56— 1.67x 10^ cm. 
The disc in KRAur is found to be untruncated or only 
marginally truncated. The inner disc radius is found to in- 
crease with the mass transfer rate, but it is always consid- 
era bly smaller than in th e case of quiescent dwarf nova (see 
e.g. ISchreiber et al.ll2003l ), when the disc is not in the steady 
state, i.e. the mass transfer rate at the inner disc radius is 
lower than the mass transfer rate from the secondary (i.e. at 
the outer disc). According to the DIM, the high mass accre- 
tion rate during dwarf nova outburst pushes the inner disc 
radius closer to the primary. The ionized steady state discs 
of nova like systems posses higher mass transfer rates than 
dwarf novae discs. Therefore, smaller inner disc radii are ex- 
pected. The maximum inner disc radius allowed by our sim- 
ulations (1.56— 1.67 X 10^ cm) are considerably smaller than 
the value 2.00 xlO^cm used f or the dwarf nova SS Cyg in 
quiescence bv ISchreiber et all l|2003l ). This is in agreement 
with the standard accretion disc theory. In summary, our 
analysis has demonstrated that a simple statistical model of 
disc turbulence is able to reproduce the observed flickering 
characteristics of KR Aur with a set of disc parameters con- 
sistent with DIM theory. The turbulence scenario is then 
a plausible mechanism and no other processes seem to be 
required to explain the observed flickering activity. 

The case of UUAqr is totally different to that of 
KRAur. The observed cut off frequency lead to a values 
lower than 0.01 in the entire range of inner disc radii and 
for all mass transfer rates. This is unacceptable even for 
discs in the low state, such as dwarf novae in quiescence. 
Concerning the power law, our simulations yield to a values 
larger than the upper limit 1.0 which is not acceptable by 
basic accretion disc theory. Therefore, no suitable disc tur- 
bulence model can succesfully reproduce the observed PDS 
of UUAqr. 

To summarize, we have two similar nova like systems, 
but the analysis of their flickering properties with our tur- 
bulence disc model lead to two completely different results. 
Is the turbulence scenario w rong for the case of UUAqr? 
iBaptista &: Bortolettd (|2008l ) found spiral structures in the 
eclipse maps of UU Aqr and they claim that the flickering is 
caused by turbulence generated by the inter action of the disc 
with the spiral density wave. The model of iDobrotka et al] 
l|2010t ) that we use in this work is based on a disc profile that 
does not take into account the presence of spiral structures. 
Furthermore, the model assumes that the largest dimension 
scale of the turbulent elements is set by the disc scale height. 
Any other structure like spiral arms can considerably affect 
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Figure 4. Evolution of the power law index (right lines) and cut 
off frequency (left lines) for KRAur as a function of mass transfer 
rate marked in units of 10^^ gs^^. 

the input distribution function of turbulence dimension. We 
propose that the presence of spiral structures generate tur- 
bulence with a totally different distribution function of di- 
mension than in a standard disc and hence our model is 
not able to describe the situation. Therefore, the turbulence 
scenario cannot be tested in this system. 

Are KRAur and UUAqr so different that one has spi- 
ral structures in the disc while the other has not ? No eclipse 
maps of KRAur are available for comparison but lKato et al.l 
(|2002| ) did not find any evidence for superhumps. Super- 
humps are brightness variations modulated with a time scale 
close to the orbital period which are caused by disc preces- 
sion. The physical phenomenon responsible for disc preces- 
sion and superhump activity is the tidal force of the sec- 
ondary star. Superhumps are triggered when the disc outer 
radius exceeds the 3:1 resonanc e radius. This is satisfied 
for mass ratios q < .35 (see e.g. IWhitehurst fc Kinel[l99ll . 
IPatterson et al.ll2005i ) . Tidal forces do not only generate pre- 
cessing discs and superhum ps but also spiral shocks in the 
disc (see e.g. ISteeghsll200ll ). According to the reported bi- 
nary parameters UUAur, with q = 0.3, is likely to be af- 
fected by the tidal influence of the secondary but not KR Aur 
with q — 0.6. Then different disc behavior may be expected 
in the two binaries. 



6.3 Distribution of flickering in the disc 

The parameter k{r) provides an estimate of the number of 
events between two adjacent disc annuli and hence it is a 
prime indicator of flickering. The lower panel in Fig.[S]shows 
that k{r) rises very steeply in the inner disc and decreases 
with increasing r. This suggests that the bulk of the events 
responsible for the flickering activity are located in the in- 
ner disc regions which is in agreement with the observations. 
iDobrotka etal] (|2010l ) explained this invoking the maximum 
dimension scale of the events. The inner disc contains very 
small turbulent eddies which are able to transfer little an- 
gular momentum. Therefore, a large number of eddies is 
needed. Conversely, the outer disc regions contain very large 
eddies which are able to transport large amounts of angu- 
lar momentum. Therefore, a smaller number of events is 
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Figure 5. Disc profile in KRAur. Radial distribution of angular 
momentum (upper panel), gradient of angular momentum (mid- 
dle panel) and correlation coefficient k{r) (bottom panel), as com- 
puted from equation [T] The very steep variations at low r values 
are due to the proximity of the white dwarf. Our simulations are 
computed at larger radii after the steep variations. 

sufficient. This is best seen when inspecting the gradient 
of angular momentum rather than the angular momentum 
itself. The top panel in Fig. [5] clearly shows that the angu- 
lar momentum increases toward the outer disc. Flickering 
is caused by the transport of angular momentum L{r) and 
hence the difference between adjacent radii is what matters. 
The gradient of L{r) decreases toward the outer disc, as 
is illustrated in the middle panel of Fig. [S] Therefore, the 
amount of L that needs to be transported decreases toward 
the outer disc and hence the quantity of turbulent eddies 
must also decrease (see lower panel in Fig. [5]). In summary, 
flickering becomes more important at low r values or the 
innermost part of the disc, as observed (see Sec [TJ . 



7 SUMMARY 

In this paper we have analyzed the flickering activity of 
two nova like systems KRAur and UUAqr. We applied 
a statistical model of flickering simulations developed by 
iDobrotka et al.l (|2010l ) based on a simple idea of angular 
momentum transport in accretion discs. The use of a steady- 
state disc model requires the disc to be in the hot ionized 
state. The two selected nova like systems are then ade- 
quate targets for this analysis. We successfully fitted the 
observed power density spectrum of KRAur using an inner 
disc truncation radius of 0.88 — 1.67 x 10^ cm and param- 
eter a ~ 0.10 — 0.40. Unfortunately, our modeling fails to 
reproduce the power density spectrum of UUAqr, proba- 
bly because of the pre s ence o f spiral structures detected by 
iBaptista fc Bortolettol (|2008l ). Such structures are expected 
in the case of UUAqr because of its low mass ratio g = 0.3. 
On the other hand, KRAur has a larger mass ratio of 0.6 
and the disc does not reach the 3:1 resonance radius where 
tidal forces from the secondary star generate spiral shocks. 
Therefore, a different disc behavior is naturally expected in 
the two systems. 

Our simulations show that the number of fiickering 
events increases at lower radii and hence concentrate in the 



inner disc regions. This behavior is explained by the gradi- 
ent of angular momentum which also increases toward the 
center of the disc. Furthermore, the duration of the turbu- 
lence events depends on the radial viscous velocity which is 
a function of the mass transfer rate and the a parameter. 
All these parameters are correlated and, although our mod- 
elling provides useful constraints, they cannot be uniquely 
determined. 



ACKNOWLEDGMENT 

This work is supported in part by the Grand-in-Aid 
for the global COE programs on "The Next Generation 
of Physics, spun from Diversity and Emergence" from 
MEXT. AD was supported also by the Slovak Grant 
Agency, grant VEGA-1/0520/10. Partly funded by the 
Spanish MEG under the Gonsolider-Ingenio 2010 Pro- 
g ram grant GSD2006-00070: "First science with the GTC" 
( http: //www.iac.es/consolider-ingenio-gtc/ I . 

REFERENCES 

Baptista R., Bortoletto A., 2004, AJ , 128, 411 
Baptista R., Bortoletto A., 2008, ApJ , 676, 1240 
Baptista R., HasweU C. A., Thomas G., 2002, MNRAS , 
334, 198 

Baptista R., Steiner J. E., Gieslinski D., 1994, ApJ , 433, 
332 

Bruch A., 1992, A&A , 266, 237 
Bruch A., 1996, A&A , 312, 97 
Bruch A., 2000, A&A , 359, 998 

Dobrotka A., Hric L., Gasares J., Shahbaz T., Martmez- 
Pais I. G., Munoz-Darias T., 2010, MNRAS , 402, 2567 

Frank J., King A., Raine D., 1992, Cambridge Astrophysics 
Series, 21 

Hameury J.-M., Menou K., Dubus G., Lasota J. -P., Hure 
J.-M., 1998, MNRAS , 298, 1048 

Kato S., Fukue J., Mineshige S., eds, 1998, Black-hole ac- 
cretion disks 

Kato T., Ishioka R., Uemura M., 2002, PASJ , 54, 1033 
King A. R., Pringle J. E., Livio M., 2007, MNRAS , 376, 
1740 

Lasota J., 2001, New Astron. Rev. , 45, 449 

Meyer F., Meyer-Hofmeister E., 1994, A&A , 288, 175 

Nauenberg M., 1972, ApJ , 175, 417 

Paczyiiski B., 1971, Acta Astron. , 21, 417 

Patterson J., Kemp J., Harvey D. A., Fried R. E., Rea R., 
Monard B., Cook L. M., Skillman D. R., Vanmunster T., 
Bolt G., Armstrong E., McGormick J., Krajci T., Jensen 
L., Gunn J., Butterworth N., 2005, PASP , 117, 1204 

Scargle J. D., 1982, ApJ , 263, 835 

Schreiber M. R., Hameury J., Lasota J., 2003, A&A , 410, 
239 

Schreiber M. R., Hameury J., Lasota J., 2004, A&A , 427, 
621 

Shafter A. W., 1983, ApJ , 267, 222 
Shakura N. I., Sunyaev R. A., 1973, A&A , 24, 337 
Steeghs D., 2001, in H. M. J. Boffln, D. Steeghs, & 
J. Guypers ed., Astrotomography, Indirect Imaging Meth- 
ods in Observational Astronomy Vol. 573 of Lecture Notes 
in Physics, Berlin Springer Verlag, Spiral Waves in Accre- 
tion Discs - Observations, pp 45 — h 



Warner B., 1995, Cambridge Astrophysics Series, 28 
Warner B., Nather R. E., 1971, MNRAS , 152, 219 
Whitehurst R., King A., 1991, MNRAS , 249, 25 
Yonehara A., Mineshige S., Welsh W. F., 1997, ApJ , 486, 
388 

Zamanov R. K., Bruch A., 1998, A&A , 338, 988 



Flickering of KRAur and UUAqr 



